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ABSTRACT 

The complex relationship between the galaxy density field and the underlying 
matter field limits our ability to extract cosmological constraints from galaxy 
redshift surveys. Our approach is to use halos rather than galaxies to trace 
the underlying mass distribution. We identify Fingers-of-God (FOGs) and 
replace multiple galaxies in each FOG with a single halo object. This removes 
the nonlinear contributions of satellite galaxies, the one-halo term. We test 
our method on a large set of high-fidelity mock SDSS Luminous Red Galaxy 
(LRG) catalogs. We find that the aggressive FOG compression algorithm 
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adopted in the LRG P{k) analysis of iTegmark et al.l (120061 ) leads to a ~ 10% 
correction to the underlying matter power spectrum aX k = 0.1 h Mpc^^ and 
~ 40% correction at = 0.2 h Mpc~^, thereby compromising the cosmological 
constraints. In contrast, the power spectrum of our reconstructed halo density 
field deviates from the underlying matter power spectrum at the < 1% level for 
k < 0.1 h Mpc~^ and < 4% at k = 0.2 h Mpc~^. The reconstructed halo density 
field also removes the bias in the measurement of the redshift space distortion 
parameter /3 induced by the FOG smearing of the linear redshift space distortions. 



Subject headings: cosmology: observations — galaxies: halos — galaxies: statis- 
tics — galaxies:elliptical and lenticular, cD 



1. Introduction 



Galaxy power spectra are an essential element in testing our cosmological models and 
studying the initial conditions of our universe. As surveys continue to grow in size, the 
statistical error bars on the galaxy power spectra continue to shrink. However, our ability 
to constrain cosmological models is limited by the systematics of relating the galaxy distri- 
bution to the underlying r natter field. For examp le, through a comparison of galaxy samples 
from 2dFGRS and SPSS. ISanchez fc Cold (120071 ) find that differing scale-dependent biasing 



between red and blue galaxies cause the cosmological constraints derived from the samples 
to differ at the ~ 2a level. These systematic problems are even more s evere in the SPSS Lu - 
minous Red Galaxy (LRG) sample. With its large effective volume ( ITegmark et al.ll2006l ). 
this sample potentially provides the tightest constraints on the shape of the linear power 
spectrum. However, because of the low number density of the LRGs and their high bias, the 
LRG power spectrum potentially requires a large nonlinear correction. As a consequence, 
several authors have encount ered problems when cor nbining the LRG pow er spectrum re- 
sults of ITegmark et al.l (120061 ) with CMB data sets. iPunkley et al.l (120081 ) find that when 
combining with the 5 year WMAP results, the best fit flrn vari es systematically with the 
maximum k from the LRG P{k) included in the analysis. When IVerde fc Peirid (120081 ) re- 
construct the primordial power spectrum with minimal assumptions about its shape, they 
find that the LRG sample has less statistical power than the SPSS MAIN or 2dFGRS, 
despite its larger effective volume. Since the authors restrict themselves to < 0.1, they 
are unable to constrain both the nonlinear correction and the primordial power spectrum. 
The difficulties in both analyses stem from the degeneracy between the power spectrum 
shape and the potentially large nonlinear correction amplitude. Fig. [1] shows that the best 
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fit iTegmark et al.l (120061 ) nonlinear correction is above the statistical error in each fc-bin for 
k > 0.07, while k < 0.1—0.15 is typically used in cosmological parameter analyses. Moreover, 
the /c-dependence of this large nonlinear correction must be accurate to avoid introducing 
systematic errors. The solution we propose in this paper is to first estimate a halo density 
field from the LRG density field. The resulting field is nearly linearly biased with respect to 
the dark matter for k < 0.2. The dashed curve in Fig. [1] shows that the nonlinear correction 
betwe en the dark ma t ter an d linear spectrum is much smaller than the correction applied 
in the iTegmark et al.l (120061 ) LRG analysis; it is also much easier to calibrate using A^-body 
simulations. 



The Halo Occupa tion Distribution (HOD) model (ISeliakll2000l : IPeacock fc Smithll2000 



Cooray fc Shethll2002l ) allows us to describe the link between galaxies and dark matter, and 
provides a framework in which to analyze non-linearities contributing to the power spec- 
trum. Pairs of galaxies can be separated as one- or two -hal o, i.e., occupying the sa me or 
distinct dark matter halos. As discussed in ISeljakl (120001 ) and ISchulz fc White! (120061 ) . scale 
dependence of the bias between galaxies and dark matter arises primarily because the one- 
halo galaxy term i s enhanced relativ e to the dark matter more than the two-halo term. To 
combat this effect, iHuff et al.l (120071 ) introduce a configuration s pace band-power estim ator 
which confines the one-halo contribution to ~ 2 — 3 h^^ Mpc. iHamann et al.l (120081 ) find 
that for measured galaxy power spectra (SDSS MAIN, 2dFGRS, and SDSS LRG samples), 
the nonlinearity can be sufficiently modeled as excess shot noise. 



In extracting the underlying linear matter power spectrum from large galaxy redshift 
surveys, there are three main approaches to disentangling the effects of nonlinear evolution in 
the matter field, nonlinear redshift space distortions, and nonlinear galaxy biasing. One can 
model the nonlinearities with a phenomenological fitting formula and marginalize over its 
free parameters when fit ting cosmological par ameters; this technique has been emplo yed in 
analyse s of SDSS LRGs fregmark et al1 l2006h as well as in 2dFGRS JCole et allboosl ). The 
work of IYoo et al.l (120081 ) instead starts with the HOD to predict the nonlinear modifications 
to large-scale clustering statistics. In principal, cosmological and HOD parameters can be fit 
simultaneously. In this forward approach, the measured large-scale clustering is compared 
with predicted clustering given the HOD parameters. The method put forth in this work is 
an inverse method: we make use of phase information in the galaxy density field to estimate 
a halo density field. The power spectrum of this field is nearly linearly biased with respect to 
the underlying nonlinear matter density field. A future exten sion of this work is to a lso use 
the phase information to reconstruct the linear density field ( lEisenstein et al.l l2007al ) . This 
combined approach would "undo" two sources of nonlinearity: galaxy bias and nonlinear 
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dynamics. 



In this paper we introduce our halo density field reconstruction method and use A^- 
body simulations to test the method. We first review previous app roaches to the analysis 
of the LRG P{k) in § H In § [3] and S we apply the CiC method of lReid fc Spergell J2008h 
to produce high fidelity mock catalogs from 42 A^-body simulations for three LRG redshift 
subsamples: NEAR, MID, and FAR. § O describes our method for reconstructing a halo 
density field. In §[6] we present the nonlinear dark matter and mock catalog power spectra and 
covariance matrices. We an alyze in detail the effec ts of the Finger-of-God (FOG) compression 
algorithm employed in the iTegmark et al.l (120061 ) analysis of the LRG power spectrum. We 
also present r e sults for the angle-averaged redshift space spectrum, which is measured in 



Percival et al.l (120071 ). We demonstrate that the reconstructed halo density field is linearly 
biased with respect to the dark matter out to k = 0.2 at the 2% level for the MID and FAR 
subsamples and at the 4% level for the NEAR sample. For k < 0.1, the traditional regime 
used for cosmo logical p arameter estimat i on, th e discrepancy is < 1%. We also show that the 
beat-coupling model of iHamilton et al.l (l2006l ) fits both the dark matter and mock catalog 
covariance matrices. Finally we examine the structure of the redshift space distortions as a 
function of k. Throughout this pa per we adopt the cosm ological parameters recommended 
from the latest WMAP5 analysis jKomatsu et al.lboosh : (l^^, fife, ^a, n„ ag, /i) = (0.2792, 
0.0462, 0.7208, 0.960, 0.817, 0.701). 



2. Background 

2.1. Summary of Previous Analyses of the SDSS LRG Power Spectrum 



Two pow er spectrum analyses were simultaneously published for the fifth data release 
of the SDSS. iPercival et al.l (120071 ) compute the monopole (i.e., angle-averaged) redshift 
space power spectrum of t he combined sample of MAIN and LRG galaxies. The method 
extends the FKP method (IFeldman et al.lll994l ) to remove differential bias from the recov- 
ered spectrum, and optimally weights the galaxies accor ding to their expecte d bias, which 
is determined from their luminosity. This approach leads iPercival et al.l (120071 ) to use a dif- 
ferent selection function for the LRGs than adopted for our analysis. We do not model the 
luminosities of the galaxies, and so precise evaluation of their method is not possible here. 
However, we do compute the monop ole redshift power sp ectrum of our mock catalogs. This 
should compare most closely to the IPercival et al.l (120071 ) method, though the weighting of 
the galaxies according to their luminosity will certainly alter the relative amplitudes of the 
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Fig. 1. — Poi nts show the error bar s divided by the bandpower for the LRG power spectrum 
pubhshed in iTegmark et al.l (120061 ) . The sohd curve shows (1 + QKhk"^) / i)- + 1-4/c) — 1, 
the fractional amphtude of the r ionhnear Qnt. correct ion to the hnear power spectrum for 
Qnl = 31, the best fit value in iTegmark et al.l (120061 ). For comparison, the dashed curve 
shows the smooth nonlinear correction fit to the dark matter power spectrum of Fig. [5l This 
is the size of the correction for our reconstructed halo density field. 
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one- and two-halo contributions to the power spectrum. 



Tegmark et aLl (120061 ) focus their analysis entirely on the LRGs. They examine three 
redshift subsamples: NEAR (0.155 < z < 0.300), MID (0.300 < z < 0.38 0), and FAR 
(0.380 < z < 0.474). We mod el each of these sampl es separately in our analysis. iTegmark et al, 
( 120061 ) use the PKL method (jTegmark et al.ll2004l ) to estimate the real space galaxy-galaxy, 
galaxy-velocity, and velocity- yelocity power spectr a. In linear theory the three spectra are 
related by (Eqns. 2 and 3 in iTegmark et al.ll2006l ): 



P'Pggik). 



(1) 
(2) 



f3 relates the amplitude of galaxy and velocity field fluctuations in linear theory (Eqn. [35]) and 
Tgy is the dimensionless correlation coefficient between the velocity field and gravitational 
acceleration field. The final estimate of the real space LRG power spectrum is a linear 
combination of all three, though it is dominated by the scaled redshift space monopole power 
spectrum. Their method assumes that the parameters relating the galaxy and velocity power 
spectra, P and Vg^, are scale independent, and they only use data with k < 0.09 h Mpc^^ to 
estimate these factor s. To make this assumption more accurate for the LRG density field, 
Tegmark et aLl (|2006f) compress FO Gs before estimating the power spectrum. The algorithm 



is detailed in ITegmark et al.l (l2004j ). but to our knowledge has not been extensively tested 
on accurate LRG mock catalogs. The alg orithm finds groups of LRGs by a Friends-of- 
Friends (FoF) algorithm (IFrenk et al.lll988l ) and then isotropizes them by scaling all radial 
separations from the group center. Two objects are considered friends when 
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-7m(l + 6c) 



-,max 



-1/3 



(3) 



with 4 = 200 and r^^^nax = 5 h'^ Mpc; the RHS of Eqn. |3] is 2.3 h'^ Mpc for the SDSS 
LRG NEAR and MID subsample number densities, and 2.8 h^^ Mpc for the FAR sample. 
For comparison, the virial radius of a IO^'^H'^Mq halo is ~ 1 h^^ Mpc (~ 2 for a 10^^h~^MQ 
halo). Since the line-of-sight (LOS) distance is compressed in the LHS of Eqn. |3l the al- 
gorithm will compress objects well beyond the virial r adius of the typ i cal ho st halo. After 
examination of the NEAR, MID, and FAR subsamples, ITegmark et al.l (120061 ) conclude that 
the expected decrease in bias and increase in clustering amplitude with redshift cancel within 
the observational errors on the subsample spectra. They therefore apply no corrections for 
evolution of the galaxy population with redshift, and compute a single power spectrum for 
the entire sample. 
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2.2. Nonlinear Power Spectrum Models 



The Qnl model was introduced in ICole et al.l (120051 ) as a fitting formula relating the 
linear and galaxy power spectra: 



Pgal{k) 



1 + QiVLfc^ 

l + Ak 



PUk). 



(4) 



Both from examining mock catalogs and the HOD model analytically, I Cole et al.l (120051 ) fix 
A = 1.4 (redshift space) and A = 1.7 (real space) to model the 'previrialization' suppression 
of power on large scales (ILokas et al.l 119961 ) . Qnl considered a nuisance parameter whose 
amplitude is set both by nonlinearities in the matter field and by the scale-dependent bias 
of the galaxy sample. 



Eisenstein et al.l (|2007bl ) examine in detail the damping of BAO features in the power 
spectrum through differential motion of pairs of tracers separated by the BAO scale. For 
our real space matter power spectrum we adopt their model 



PsmeUk) = Pun{k)e-'"/''BAO + 



no wiggles 



{k)(l 



(5) 



-Pno wiggles (fc), defined by Eqn. 29 of lEisenstein &: Hul (1l998l ). is a smooth version of Punik) 
with the baryon oscillations removed. This means the baryon oscillations remaining in 



Psmear{k) are damped completely for k ^ ksAo and unaltered for k <^ ksAo- In § 16.11 we 
find kBAo = 

0.14 h Mpc"^ for our real spa c e matt er power spectrum, in agreement with 
the value reported in §6 of lEisenstein et al.l (l2007bl ): they also find that Eqn. [5] is a good 
approximation in redshift space, but with smaller ksAo than in real space. 



Percival et al.l (120071 ) assume that both P{k) and galaxy bias are linear, and show that 
the resuhing fits to fi^ for 0.01 < k < 0.06 h Mpc"^ and 0.01 < k < 0.15 h Mpc"^ 
are discrepant at the 2-3a level, demonstrating the need for nonlinear modeling to extract 
robust cosmological inf ormation from the obs erved galaxy power spectrum. The LRG power 
spectrum measured in iTegmark et al.l (120061 ) requires a large nonlinear correction (see our 
Fig. [1]). They adopt the model 

1 + QNLk^ 



PLRcik) = b Psmearik, ksAO = 0.1) 



(6) 



l + Ak 

with A fixed at 1.4 (M. Tegmark, private communication). Their best fit value for the nui- 
sance parameter is Qnl = 30. This model is only tested on power spectra of mock LRGs 
in real space and with no satellite LRGs. We will compare this model to our more realistic 
mock catalogs in redshift space. 
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2.3. Simplest HOD Model for the LRG Power Spectrum 



As we describe in more detail in § 14. ![ the HOD model sp ecifies P(Nr,RG\M), th e 
probability that Nlrg LRGs occupy a halo of mass M. Following Cooray fc ShethI (120021 ) . 
the nonlinear power spectrum can be separated into the contribution from pairs of galaxies 
occupying the same dark matter halo, P^^{k), and p airs occupying di s tinct dark matter 
halos, P^^ik). In the large scale limit of Eqn. 128 in ICooray fc ShethI (120021 ) . where the 
effect of finite size of halos can be safely neglected and the bias with respect to the dark 
matter can be approximated as linear, the LRG power spectrum is approximately 



PLRcik) 



Ih 

LRG 



PhRGi^) + PlRgO^) 



Plrg^^) 



^LRG^DMik) 



''^LRG 



(7) 



(9) 



where n{M) is the halo mass function, Blrg is the large scale bias between the galaxies and 
matter, and PoAiik) is the nonlinear dark matter power spectrum. This simple model as- 
sumes that the effects of satellite galaxies are simply to reweight halo pairs traced by central 
galaxies in P^'^, and to add a simple shot noise term given by Eqn. [81 As we show below, 
this division of pairs works well in real space, and PhRGik) — b^RcPDuik) is well described 
by P^^. However, pairs of galaxies in redshift space cannot be cleanly described by this 
separation. For our HOD parameters, the halo velocity dispersion causes satellite galaxies 
to shift by (a^)^^^ ~ 9 h^-^ Mpc (with large tails) along the LOS. This shifting causes power 
to be shuffled between scales and eve n the largest scale modes along the LOS to be damped 
by the FOG smearing; iTinkeii (120071 ) provide a detailed treatment. Our approach will be to 
circumvent many of the difficulties in redshift space by eliminating the FOGs due to satellite 
galaxies from the density field before computing PLRG^k). 



2.4. Modeling the Covariance Matrix 

The covariance matrix for a set of band powers -P(A;j=i..Ar) is defined by 



Cij 



[{P{h) - P{h)){P{k,) - P{k^)) . 



(10) 



Hamilton et al.l (l2006l ) showed that the largest nonlinear contribution to the covariance ma- 



trix is a beat-coupling term proportional to the power on the largest scale of the survey. 
In any real survey of finite volume, the observed Fourier amplitudes are convolved with the 
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survey window function Ws'. 



dk' 

(2^ 



The beat coupling contribution arises because neighboring Fourier modes (5(k)5(— k — e) are 
coupled by nonlinear growth to the beat mode S{e). When the DC mode of the survey is 
positive, all modes are amplified; when it is negative, all modes are suppressed. This terra 
can be large since the linear power spectrum drops so sharply with k. [Hamilton et al.l (120061 ) 
emphasize that this term does not contribute to the covariance of power measured from 
ensembles of traditional periodic box simulations, where the band power ki is averaged over 
Ni complex modes by 



(12) 



Here, the beat mode is k,;j — kjj 



0, the DC mode. However, since our simulations allow 
the DC mode to vary (jSirkd l2005l ). we will capture this term. We therefore model our 
covariance matrix as the sum of the usual diagonal Gaussian and shot noise contributions 
and the beat-coupling contribution: 



ij,model 



1 



+ ^) 5J + A(3,,,,R^P{k:)P{k^)5lc{z) 



n 



(13) 



Here 5^, is the Kronecker delta, and S'j^ciz) 



PL{k = 0,z = 0)/Vs.my<D\z)/D' 



= 0) is 

the variance of the DC mode linearly evolved to re dshift z in the simulat ion volume Vsim = 
(1Sirkdl2005l ). In perturbation theory Ra ~ 2.62 (jHamilton et al.ll2006l ). We introduce ashot 
and Pbeat to allow for excess shot noise and variation in the amplitude of the beat coupling 
term, though we expect both parameters to be ~ 1. We show in § E]that Eqn. [T^ provides 
a good model for both the dark matter and LRG covariance matrices. For the dark matter, 
the shot noise contribution is negligible. 



Neyrinck et al.l (120061 ) showed that Poisson fluctuations about the mean halo mass func- 



tion introduce variance in the amplitude of the one-halo contribution to the dark matter 
power spectrum that can dominate the covariance matrix in the nonlinear regime. Since 
our reconstruction of the halo density fleld seeks to eliminate the one-halo contribution from 
galaxies, we expect the covariance to be smaller for the reconstructed halo density field than 
for the original galaxy sample. 
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2.5. Error Estimates on P{ki) and Cij 

For the dark matter, errors on the bandpowers are straightforward. They are simply 
the diagonal terms of the inverse covariance matrix divided by the number of simulations. 
We use the model covariance matrix in Eqn. [13] to compute the inverse: 

^P{h) = 1^7 i^model)ii- (14) 

Of course, the beat coupling leads to off-diagonal terms in both the covariance matrix and 
its inverse; these terms must be included when estimating errors on model parameters. 



We also estimate an error on our estimate of the covariance matrix. We ignore correla- 
tions between the Cj/s in our error estimates, though they are certainly present. 



1 



Cij AT . 



'{{p{h) - p{mp{k,) - p{k,)) - . (15) 



The situation for the mock catalogs is more tricky. For each of the s = 1, .., Ngim iV-body 
simulations we produce m = 1, .., Nmocks mock catalogs using our fixed HOD parameters to 
reduce the shot noise contribution. We define 

Cij,HOD = {{Psmih) - Ps{hd){Psm{k,) - Ps{k,))) (16) 
C,,,red = {{Psih) - P{h)){P,{kj) - P{k,))) (17) 

Qj,tot = {{Psmih) - P{k)){PsUkj) - P{k,))) (18) 

HOD + Cij,red- (19) 

Here Psm{ki) denotes the band power for mock catalog m populating simulation s, Ps{ki) 
denotes the band power in a single simulation s averaged over the m mock catalogs populating 
s, and P{ki) denotes a band power averaged over the entire set of Nmocks x Ngim catalogs. 
Cij, tot is the covariance matrix for this set of mock catalogs. Cij^HOD is the variance introduced 
by sampling the same matter density field with different mock catalog realizations and Cij^red 
is the reduced covariance of the power spectra from each simulation after averaging over 
Nmocks catalogs in each simulation. The expected error on bandpower P{ki) is then 



P{ki) = {{Psm{h)) (20) 

1 

nI 



^P(ki) — model, red) a- (21) 



We again neglect the covariance of the covariance matrix elements to estimate the errors on 
Cij^tot, Cij^red, and Cij^HOD from their variance over the simulations. These errors are only 
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used for fitting the parameters ashot and (3beat, so inaccuracies in the errors should make our 
estimates noisier, but not biased. 



Simulations 



We use the publicly-available Tree-Particle- Mesh (TPM) code ( iBode fc Ostriker 120031) 
to run 42 periodic box simulations with the parameters selected in Chapter 3 of iReid l2008h : 
Lbox = 558 Mpc and Np = 512^, corresponding to a particle mass Mp = 1.43 x lO^^Af©. 
This resolution ensures that there are at least 50 particles per halo populated with an LRG 
for our best fit HOD parameters. We allow the slight change in the halo mass function in 
our lowest mass bin to be absorbed by a slight change in the best fit HOD parameters for 
central galaxies. As outlined in Appendix A, we performed a battery of tests to ensure that 
the large and small scale clustering and velocity statistics were unaffected by our choice of 
mass resolution for the simulation. 



Initial conditions for the simulations were generated with the publicly available ic code 
(jSirkdl2005r). We updated the code to use the mtl9937 "Mersenne Twister" random number 
generator (IMatsumoto fc Nishimuralll998l ) to select the initial density mode realizations, ic 
was designed to generate initial conditions for an ensemble of periodic box simulations that 
will match real-space statistical properties such as the mass variance in spheres and ^(r). 
For a periodic box simulation with comoving side length Luni at z = oo, the initial modes are 
drawn from a convolved power spectrum Pi^^iik) for which PLu„i{^) / ^tni variance of 

the DC mode in a volume equal to the initial periodic box volume. The DC mode is assumed 
to evolve linearly, and so its effects can be mimicked to first order in Lagrangian perturbation 
theory by a slight change in the cosmological parameter s of th e box a s a function of the value 
of the DC mode in the realization (see Eqns. 19 - 22 of ISirkd (120051 ) ). The relation between 
the scale factor in the box cosmology, abox, and the scale factor in the universe cosmology, 

O^ii.rt.i^ IS 



1 Dun 

" 3~d: 



-A, 



(22) 



where Aq is the amplitude of the DC mode at auni = 1, and the growth function Duni is 
evaluated using the universe cosmology. In this scenario, each simulation in the ensemble 
represents an equal initial comoving volume, L^^nv Once the universe has evolved to scale 
factor auni, the comoving size of the box is Lbox = Luni^cibox / CLuni) ■ Because overdense regions 
expand more slowly than underdense regions, we must weight the simulations by {abox / duni)^ 
to obtain volume-averaged quantities. This improvement over a set of simulations for which 



the DC mode is precisely is crucial to the analysis in this paper, since we wish to extract 
accurate matter and halo power spectra as the density field enters the nonlinear regime, 
along with accurate estimates of band power covariances. 



We output the dark matter particle positions and velocities at time steps nearest to 
universe redshifts znear = 0.235, zmid = 0.342, and zpa r = 0.421, the galaxy - weighted 
mean redshifts of the three redshift subsamples analyzed in lTegmark et al.l (120061 ). Eqn. [221 
shows that for a snapshot of the universe at fixed redshift, each simulation should be output 
at a different Zbox that increases with Aq. Since the TPM code outputs particle data at the 
nearest whole time step to the desired redshift, the difference between the target and output 
redshift will vary slightly with A^. We correct for this effect exactly in the linear regime by 
scahng the power spectra from each simulation by {D {ztar get) /D{z output))"^ before comput- 
ing averaged quantities. Without this scaling, we get the same average dark matter power 
spectra; changes in the covariance estimates are well below 1%. Any nonlinear corrections 
to this scaling will have an even smaller effect, and so we safely neglect them. 



3.1. Calculating Power Spectra 



Because of the low number density of LRGs, the shot noise correction is large compared 
to the P{k) of the continuous density field sampled by the LRGs. We use an FFT with 
= 512^ points, and present the power spectrum out to k = 0.4 h Mpc~^ = O-lAk^, 
where kj^j = n/Agrid is the Nyquist frequency. Comparis on with an iV = 1024^ grid s h owed 
agreem ent at the ~ 10~^ level. Follow ing the work of iHockney fc Eastwood! (jl98ll ) , iJing 



(120051 ) . and lBaugh &: Efstathioul (119941 ). mock LRGs are distributed on the FFT grid using 
the t riangular-shaped cloud (TSC), and the power spectrum is estimated at each k by (jjing 
2005h 

\SlRG,FFt(^)\'^ — Pshot,TSc(^) 



P(k) 



shot,TSC[ 



W' 



TSC{ 



n 



1 — sin^ 



Tckj 



smc 



irki 



2k 



N 



2kN 

3 



15 



sm 



nki 



2k 



N 



(23) 
(24) 

(25) 



where Slrg,fftO^) is the FFT of the LRG overdensity field and ki are the Cartesian com- 
ponents of k. Eqn. [23] corrects for FFT aliasing in the limit that P(k + kN) is shot noise 
dominated, which is the case here. Finally, we average the values of Eqn. [23] over k-bands 



with Akbox = 0.0113. 



3.2. Halo Catalogs 



We use the spherical overdensity (SO) halo finder code described in [Tinker et al.l (120081 ) 
with A = 200p5 to generate halo catalogs. Fig. [2] shows the rat io of the volume-we ighted 
SO mass function in our simulation set to the fitting function of [Tinker et al.l (120081 ) down 
to M = 7.15 X IO^^Mq. There is good agreement at the level of accuracy claimed for their 
analytic fits (~ 5%). The overabundance of halos in the lowest mass bins is likely due to 
the small number of particles per halo in those bins. This slight modification of the mass 
function is unimportant for our purposes, since there are no satellite LRGs in these mass 
bins, and changes in the mass function are degenerate with changes to the probability of 
hosting a central galaxy as a function of halo mass, Ncen{M). 



4. Mock Catalogs 



We use the technique described in iReid fc Spergell (120081 ) to produce mock catalogs for 
the NEAR (0.155 < z < 0.300), MID (0.30 < z < 0.380), and FAR (0.380 < z < 0.474) 
redshift subsamples of iTegmark et al.l ( l2006l ). These samples contain roughly equal numbers 
of galaxies, and their power spectra have roughly the same amplitude on l arge scales. We 
first r eview our HOD modeling assumptions, described in more detail in iReid &: Spergel 
(l2008h . 



4.1. HOD Model 



The Halo Occupation Distribution (HOD) model assumes that the probability P{NLFtG\M) 
of Nlrg LRGs occupying a dark matter halo of ma ss M at redshift z depends only on the 
halo mass (for a review, see ICooray &: Shethll2002l ). However, the application of the HOD 
formalism requires that we make several further assumptions about P{Nlrg\^)- Detailed 
studies of dark r natter halos and subh alos suggest a division of galaxies into central and 
satellite galaxies (IKravtsov et al.l 120041 ). The central galaxies are assumed to sit at the halo 
center, consistent with the observation that most (~ 80%) of the brightest cluster LRGs are 
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Fig. 2. — The rati o of the volume- weig hted SO mass function from our simulations with the 
spline fit given in [Tinker et al.l (120081 ) for the NEAR (solid line), MID (dashed line), and 
FAR (dotted line) in bins of AlogioM = 0.1. The agreement is within the stated accuracy 
of the iTinker et al.l toO^ spline fit (~ 5%). 
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found within 0.2rnr of the center of the cluster potential well as traced by X-rays. Satellite 
galaxies occur in the more massive halos already containing a central galax y. In high resolu- 



tion s imulations they can be directly associated with dark matter subhalos (IVale &: Ostriker 



20061 ): here we will assume they have the same distribution as the halo dark matter. We will 
use these functional forms with the five free parameters Mmin, criogM, Mi, Mcut, and a to 
describe the mass dependence of the average halo occupation as a function of halo mass M : 



(Ncen) = \ 



{N{M)) = (iV,,„) (1 + {N,at)) (26) 
logioM - logioMrm 



1 + erf 



C^logM 



(27) 



Halos are populated with a central galaxy with probability Ncen{M). Central galaxies are 
placed at the center of their host halos and assigned the peculiar velocity of their halos. 
Halos with a central galaxy are populated with Ngat galaxies, where P{Nsat\N{M)) is drawn 
from a Poisson distribution. The position and velocity of a satellite galaxy is taken to be 
that of a randomly selected dark matter particle halo member. We assign comoving redshift 
space position s to an object in our mock catalogs using the conversion at Zbox, the redshift 
at which the simulation data were output: 

S = XlOS + (1 + Zbox)Vp/H{Zbox) (29) 

where xlos is the comoving distance along the line of sight in real space. 



4.2. Measurements of Ncic from SDSS and HOD parameter results 



The HOD mod el specified in ^14. II dete rmines the number of groups with Usat = 0, 1, 2, ... 
satellite galaxies. In lReid fc Spergell (120081 ) we present the measurement of the SDSS Counts- 
in-Cylinders (CiC) group multiplicity function Ncic{nsat), calibrate its relation to the true 
group multiplicity function, and use a maximum likelihood analysis to derive the parameters 
of Nsat{M) in Eqn. [28l The parameters of Ncen{M) in Eqn. [27] are constrained by the ob- 
served LRG number density and large scale clustering amplitude. To derive HOD parameters 
for the mock catalogs used in this paper, we follow the technique of iReid &: Spergell (120081 ) 
exactly, but we must first measure Ncicinsat) separately for our three redshift subsamples. 



In Table[T]we report our measurement of Ncicinsat) from the SDSS DR4+ LRG sample 
for the NEAR, MID, and FAR redshift subsamples. We introduce a few minor changes in 
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this measurement from lReid &: Spergell (120081 ) . First, we include a boundary set of galaxies 
in the redshift direction, since FOGs may widely separate nearby pairs of LRGs; we verified 
that this has a small effect on the resulting CiC group multiplicity. This buffer is Az = 0.007 
at the low redshift end of the sample and Az = 0.009 at the high redshift; that is, equal to 
the maximum redshift separation for a CiC pair of galaxies at the sample boundary. This 
causes the redshif t ranges of the subsarn ples for which we measure Ncicip-sat) to be slightly 
different than the lTegmark et al.l (120061 ) NEAR, MID, and FAR samples. Because the color 
cuts used to select the LRG sample produces a complex radial selection fun ction, the effec- 
tive n umber density of the two samples is slightly different. As discussed in iReid fc Spergel 
(120081 ) . the main difference betw 'een our sample num ber densities (Column 4 of Table [2]) 
and the result of integrating the IZehavi et al.l (120051 ) model for the redshift dependence of 
^lrg{^) (Column 5 of Table [2]) comes from our careful inclusion of objects from the imaging 
sample. In the MID sample, this is an 8% increase in the number density. We have verified 
that once the imaging galaxies are accounted for along with the difference in redshift range of 
the NEAR and FAR samples, our simple nu mber density estima te ulrg = N sample /Vsampie is 
in agreement with the expectations from the lZehavi et al.l (120051 ) model. Second, our redshift 
indicator used for objects without spectra r nust be altered, as the 4000 A break moves from 
the g to the r band at z ^ 0.4 (see Fig. 4 of lEisenstein et al.ll200ll ). For our FAR sample we 
use the r — i color as a redshift indicator. We use the observed Ncic{nsat) for the NEAR, 
MID, and FAR subsamples to derive the subsample satellite HOD p arameters Mr.ut, Mi, an d 
a in Table [2] using the max imum likelihood technique presented in iReid fc Spergell (12008! ). 
Since iTegmark et al.l (120061 ) find that the NEAR, MID, and FAR subsamples have consistent 
P(/c)'s, we vary aiogM to r natch each redshift sub sample to the large scale amplitude of the 



combined P{k) reported in lTegmark et al. 



the agreement of the large scale P(k) with 



20061). M^,r, is determined by Usim- Fig. |3] shows 



Tegmark et al.l (120061 ) when we apply their FOG 



compression algorithm to our mock catalogs. 



5. Reconstructing the Halo Density field 



In the CiC technique detailed in iReid fc Spergell (120081 ) . two galaxies are considered 
neighbors when their transverse comoving separation satisfies Ar_L < 0.8 h^^ Mpc and their 
redshifts satisfy Az/{1 + z) < Avp/c = 0.006. A cylinder should be a good approximation 
to the density contours of satellites surrounding central galaxies in redshift space, as long as 
the satellite velocity is uncorrelated with its distance from the halo center and the relative 
velocity dominates the separation of central and satellite objects along the line of sight. 
Galaxies are then grouped with their neighbors by a FoF algorithm. The reconstructed halo 



Tabl e 1. The final Nninin .^at) group multiplicity function following the method in 
Reid fc Spergdl feoosh for our NEAR, MID, and FAR LRG subsamples. 



risat 


^CiC,NEAR(n) 


^'^CiC.MIDCn) 


.'VciC,PAR("') 





22921.71 


24537.81 


19109.71 


1 


1372.63 


1301.29 


664.28 


2 


170.01 


153.40 


61.94 


3 


41.85 


25.59 


6.08 


4 


15.16 


9.07 


2.04 


5 


2.11 


2.04 


1.00 


6 


1.01 


1.01 


0.00 


7 


1.02 


0.00 


0.00 


8 


0.03 


0.00 


0.00 



Table 2. Mock catalog parameters. Masses in units of lO^^M© ^-re for SO halos with 
A = 200pfe; number densities are in units of 10"'' {h~^ Mpc)~^. zsDss,CiC is the redshift 
range of the SDSS LRG subsample and Zgim is the average simulation Zuni- ngim is the 
mock catalo g number density wh ile Umodei is the expected number density obtained by 
integrating the IZehavi et al.l (120051 ) model over the zsDSS,CiC range. P^'* is computed from 

Eqn. [8] in units of (h~^ Mpc)^. 



Sample ZSDSS.CiC ^sim risim, n,nodel O-iogM Mmin M^ut Mi a P^^ 

NEAR 0.162 - 0.300 0.2347 1.05 0.982 0.6 0.78 0.49 5.87 1.16 1759 

MID 0.300 - 0.380 0.3420 0.962 0.890 0.6 0.76 0.57 6.29 1.05 1480 

FAR 0.380 - 0.465 0.4216 0.470 0.422 0.9 2.27 1.38 5.97 0.78 2312 
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Fig. 3. — kP{k) for the NEAR (dotted curve), MID (dashed curve), and FAR (soli d curve) 
mock galaxy samples with the FOG compression algorithm of iTegmark et al.l (l2006l ) applied 
compared with their observed power spectrum (points with large error bars). Our NEAR 
and MID samples are indistinguishable, so we show error bars only for the MID and FAR 
samples. Our error bars are de rived from the d i agona l elements of the inverse Cij^red matrix. 
We scale the real space P{k) of iTegmark et al.l (120061 ) by (0.8)^^ to approximate the redshift 
space monopole (i.e., angle averaged) power spectrum; this is the relation between the real 
and redshift monopole spectra if one neglects the small contributions from the quadrupole 
and hexadecapole detailed in their §A3. For comparison, the long dashed curve shows our 
simulation initial conditions drawn from the convolved power spectrum according to the 
ic algorithm, and the smooth dot-dashed curve shows the linear power spectrum for the 
cosmological parameters adopted in this study. 
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density field is defined by the superposition of the centers of mass of the CiC groups. The 
CiC parameters were established as a balance between completeness and contamination in 
identifying pairs of galaxies at an early stage of this work based on FoF halo catalogs. While 
the parameters and details of the method described here are sufficient for approximately 
recovering the halo density field power spectrum, the method could almost certainly be im- 
proved to more accurately recover groups of galaxies residing in the same dark matter halo. 



6. Results 
6.1. Dark Matter 

6.1.1. Matter Power Spectrum 

Fig. m shows the ratio of the matter power spectrum to the power spectrum of the sim- 
ulation initial conditions scaled by the expected linear growth D'^{z) for our three redshift 
samples; the results at different redshifts are highly covariant since they are measured from 
the same set of simulations at relatively small separations in time. Error bars are computed 
from the inverse of the model covariance matrix in Eqn. [13] with /?beat = 1 fixed. 



As expected, the nonlinear correction grows as the redshift decreases. The nonlinear 
evolution generates power that smoothly in creases with fc, ari d also damps the baryon os- 
cillations. The halofit nonlinear correction (jSmith et al.ll2003l ) (dotted curve, evaluated at 
zmid) underestim ates both the baryon wiggle su ppression and the amplitude of the smooth 



Crocce fc Scoccimarrol (120081 ) also find a disagreement between their A^- 



mcrease m power; 
body simulation results and the halofit nonlinear correction. 



Our fitting function to the nonlinear matter power spectrum Pdm(^) aims to capture 
both the smearing of th e acoustic peaks and a smooth increase in power with k. We follow 
Eisenstein et al.l (l2007bl ) in defining Psmear{k) (Eqn. [5]), but adopt a slightly different fitting 
function for the smooth correction. Our model is 



Pdm ; kBAo) (^0 + o.ik + a2k'^ + as/c^) . 



(30) 



For the MID sample with 35 bandpowers and 5 parameters, = 21. W e find ksAo = 0.14, 
in good agreement with the values reported in lEisenstein et al.l (l2007bl ). In Fig. Owe show 



that the -Psmear(^) term accounts for the baryonic features in PDj\,/(fc), and the polynomial 
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in k adequately fits the smooth correction for the MID sample, while halofit underestimates 
the smooth correction by ~ 4% at = 0.2. The NEAR and FAR fits are similar. In Table [3] 
we list fits for the NEAR, MID, and FAR power spectra out to a maximum k of 0.2 and 
0.4 h Mpc^^. When k^axjit = 0.2 h Mpc~^, only the first three terms in the polynomial 
expansion are necessary for a good fit. 



6.1.2. Covariance Matrix 

Fig. E] shows diagonal elements of the normalized covariance matrix, Cu/ P{kiY . We 
estimate the errors from the diagonal variances (Eqn. [TSll : these may not capture the true 
errors since off-diagonal elements should be present in the 8-point function as well. Never- 
theless, when we use these error estimates to compute for the model in Eqn. [131 we find 
= 1600 for 1296 degrees of freedom (0 < A; < 0.4); if we restrict the covariance matrix 
to the 196 elements with both k bands between 0.056 and 0.21, we find = 262. In this 
case there are no free parameters and we deem the model a reasonable fit. If we allow the 
amplitude of the beat coupling term to vary, we find a best fit value 0.96 for the full matrix 
(x'^ = 1564) and 0.84 for the 0.056 < k < 0.21 subsample (x^ = 218). We note that for 
the DC mode realizations of our 42 simulations, the variance is a factor of 0.83 lower than 
the expected variance (in agreement with the expected random variation for a single mode, 
^sim = 15%). We conclude that Eqn. [13] is an excellent fit to our dark matter covariance 
matrix. 

Another point of interest in Fig. [S] is that the inverse of the diagonal elements of the 
inverse covariance matrix are nearly equal to the Gaussian expectation (dashed curve). This 
means that while the beat-coupling does not introduce additional errors on the measurement 
of the bandpowers, the measured values will be covariant. However, for a beat-coupling term 
in the form of Eqn. [121 only information on the overall amplitude of P{k) is lost. However, 
large scale structure analyses traditionally marginalize over the amplitude of P{k), since as 
and bias are degenerate. 
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0.1 

k h/Mpc 



Fig. 4. — The ratio of the volume-averaged matter power spectra at zfar = 0.421, 
zmid = 0.342, and znear = 0.235 to the input power spectrum scaled by D'^{z)/D^{z = 0); 
the amplitude of the nonlinear power spectrum increases as the redshift decreases. The non- 
linear correction to th e matter power spe ctrum measured from our simulations is larger than 
expected from halofit (ISmith et al.ll2003l ). shown as the dotted curve and evaluated at zmid- 
We also overlay -P/c'(^)/-Pno wiggles (^) + 0.4 (solid curve oscillating about 1.4) to indicate the 
location of the baryon wiggles in the initial power spectrum. The nonlinear evolution gener- 
ates power that smoothly increases with k and damps the baryon oscillations. Error bars are 
shown only for the MID sample for clarity and are estimated from the inverse of the model 
covariance matrix defined in Eqn. [131 where the shot noise is negligible and /?feeai = 1; the 
beat-coupling term in Eqn. [13] correlates the bandpowers. 
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Fig. 5. — PoMik) / Picik) for the MID sample (solid curve, as in Fig.H]), PnAiik) / Psmea.r{k) for 
ksAo = 0.14 (points with error bars), and our polynomial fit 1.026 — 1.199fc+ll. 06^^—8.426/;;'^ 
(dashed curve). We also overlay -P/c (^)/-Pno wiggles (/i^) + 0.4 (solid curve oscillating about 1.4) 
to indicate the location of the baryon wiggles in the initial power spectrum. Using error 
bars derived from the inverse of the model covariance matrix defined in Eqn. [131 = 21 
for 35 bandpowers and 5 parameters. Comparison of the solid and dashed curves shows that 
Psmca.r{k) adequately models the BAO features , and that our thi rd-order polynomial in k is 



sufficient for k < 0.4 h Mpc~^. The halofit ( Smith et al.ll2003l ) correction (dotted curve) 
shows better agreement with the smooth portion of the nonlinear correction, though halofit 
underestimates the nonlinear power by 4% at A; = 0.2 h Mpc~^. 
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Fig. 6. — The points with error bars show the measured diagonal elements of the covariance 
matrix. Error bars are estimated from the variance of Ca across our simulations. The dotted 
curve shows the Gaussian prediction for the diagonal covariance, the inverse of the number of 
complex modes in each linearly-spa ced fc-bin. The solid cu rve shows the sum of the Gaussian 
term and the beat-coupling term in [Hamilton et al.l (120061 ). which provides a good model for 
the full covariance matrix. The dashed curve shows the inverse of the diagonal elements of 
the inverse covariance matrix. We find that while the errors on the bandpowers are very 
close to the Gaussian expectation, the beat-coupling terms in the covariance matrix mean 
the bandpowers are correlated. 
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6.2. Mock Catalogs 

6.2.1. Power Spectra 
Fig. [3] shows the agree ment between ou r NEA R, MID, and FAR power spectra and the 



measurements presented in iTegmark et al.l (120061 ) . We used this comparison to set by eye 



the large scale bias of our mock catalogs through aiogM, the width of Ncen{,M) in Eqn. [271 
The small discrepancy between the normalizations of the FAR and other catalogs at small k 
could be eliminated with a slight variation in a\o^M,FAR- Ho wever, the FAR samp le clearly 



has a different shape, and is consistent with trends in Fig. 6 of iTegmark et al.l (120061 ). though 
in the SDSS sample the significance of these trends is less clear. Our FAR sample has ~ 3% 
more power at A; = 0.09 and ~ 10% more power at A; = 0.2 than the NEAR and MID samples. 



Having established the agreement between our mock catalogs and the ITegmark et al. 



(120061 ) observed PLRcik), we now attempt to isolate the sources of nonlinearity in our cata- 
logs and demonstrate that the power spectrum of the reconstructed halo density field is the 
best tracer of the underlying linear spectrum. In Fig. [7]we analyze in detail the MID sample, 
while Figs. E] and [HI show that the NEAR and FAR sample behave similarly. The first major 
result is that there is no detectable deviation from a constant bias between the central LRGs 
and the dark matter in real space for k < 0.1, and the deviation atA; = 0.2is<2% for the 
MID and FAR subsamples; the discrepancy for the NEAR sample is ~ 1% at /c = 0.1 and 
5% at A; = 0.2 (solid curves with the lowest amplitude as A; — 0). The lower left panel of 
Fig. El demonstrates that the main effect of using redshift space coordinates for the central 
LRGs is the further damping of the BAO signatures. There is also a nonmonotonic but 
smooth variation that is < 4% between A; = and k = 0.4. In the k < 0.2 h Mpc~^ regime, 
the direction of the deviation is opposite of that between the real space central LRGs and 
underlying dark matter, so that the redshift space central LRG power spectrum is nearly 
linearly related to the real space matter power spectrum at A; < 0.2h Mpc~^. 



The inclusion of satellite LRGs in real space has two effects. The linear bias is increased 
because the satellites upweight only the most massive, more highly biased halos traced by 
the LRGs. Secondly, the satellite galaxies add a shot noise given by Eqn. [HI The upper 
right panel of Fig. [7| shows that the difference between the P{k) including satellites and 
P{k) of the central LRGs only is well described by these two effects, with a relative bias 
brei = 1.042 and P^^ ~ 1460 {h~^ Mpc)^. This latter value is in good agreement with 1480 
{h~^ Mpc)^, the value computed directly from Eqn. [SI using our input HOD parameters. As 
k — > 0.4, the amplitude of the difference slightly diminishes due to the width of the satellite 
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LRG density profile within the halos. P^^ can be ~ 7 — 11% of PLRcik) at k = 0.1, and 
~ 15 — 23% at k = 0.2 (see Table [2]). If uncorrected for, this will be the dominant source 
of nonlinearity for LRGs owing to their low number density, which appears in the denomi- 
nator of Eqn. [HI Moreover, since the expected P^'^ varies across the NEAR, MID, and FAR 
samples, the power spectrum shape will necessarily vary with redshift, as is evident in Fig. [31 



We now consider the implications of the FOG-compression algorithm of iTegmark et al. 



(120061), which was designed to recover the real space PhRcik). Figs. [7| through [9] demonstrate 
that their compression algorithm (Eqn. [3]) is overly aggressive and adds more one-halo power 
than is present in the real space mock c atalog (dash-dot curves). This is the main source of 
the large nonlinear correction req uired inlTegmark et al.l (120061 ). However, Fig. [TUl shows that 
the nonlinear correction model of iTegmark et al. J2006h in Eqn. [6] accurately describes their 
resulting power spectrum for each of the redshift subsamples, so that the r esulting cosmo 
l ogica l constraints should not be biased. The best fit ACDM - Qnl model in lTegmark et al. 



( I2OO6I ) (long-dashed curve of Fig. [TO!) is consistent with being a weighted average over the 
NEAR, MID, and FAR subsample mock catalog power spectra. The dot-dashed curve in 
Fig. [TU] shows the ratio of the linear P{k) for their best fit cosmological parameters with 
Pjc{k) from our simulations. These two models differ at a level larger than the statistical 
error bars shown in Fig. [H even if we restrict the analysis to k < 0.1, and certainly if we 
can use the measurements between k = 0.1 and k = 0.2. However, the Qnl parameter is 
degenerate with this difference in power spectra, so that with this nonlinear correction model 
we cannot distinguish the two linear spectra. 



In redshift space the effect of the satellite galaxies on the angle-averaged power spec- 
trum is more complicated. At high k the ratio P{k)/ PijM{k) turns over as the suppression 
of power by the FOGs becomes more important than the addition of one-halo power. In 
the bottom right panel of Fig. [7] we see that the difference between the all LRG and cen- 
tral LRG power spectra is no longer well-approximated by the constant P^'^ ~ 1480 {h~^ 
Mpc)^ from Eqn. [8] The one-halo pair separation in redshift space is large (~ 9 Mpc 
but with broad tails), so that the neat division between one-halo and two-halo pairs in real 
space no longer holds, and power will be transferred between /c-bands in a complicated way. 
The reconstructed halo density field method removes the satellite galaxies b efore comput - 
ing P{k), which allows us to avoid the necessarily complicated FOG modeling (lTinkerll2007l ). 



The power spectrum of the reconstructed halo density field is very similar to the redshift 
space power spectrum of the central LRGs for the NEAR, MID, and FAR subsamples, as one 
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would expect if the reconstruction is sufficiently accurate. The slight difference is evident 
in the bottom left panel of Fig. [71 and may be because our method puts the reconstructed 
halo at the center of mass of the CiC group, which will leave residual FOG smearing. An 
improved method may be to put the reconstructed halo at the position of the brightest LRG, 
which should be closest to the halo center and have a smaller velocity with respect to its 
host halo. For < 0.1 the reconstructed halo density field does not detectably deviate from 
the shape of the underlying matter power spectrum for the NEAR and MID samples; for 
the FAR sample, the deviation is ~ 1%. In contrast, the FOG-compressed sample deviates 
from a constant bias at the 6%, 7%, and 10% level between k = 0.05 and k = 0.1. The 
reconstructed halo density field has only small deviations from a scale-independent bias out 
tok = 0.2: the deviation is 4%, 2.8%, and 2.5% for the NEAR, MID, and FAR sample. The 
systematics out to k = 0.2 should therefore be small enough for use in cosmological anal- 
yses. In constrast, the FOG-compressed spectra deviate from a constant bias at the 19%, 
20% and 30% levels for the NEAR, MID, and FAR samples; this is a factor of ~ 5 larger 
than the statistical errors on the bandpowers in iTegmark et al.l (120061 ) between k = 0.1 and 
k = 0.2. The nonlinear correction must therefore be extremely well-calibrated to extract any 
cosmological information in this regime. 



We have shown that the nonlinear correction between the reconstructed halo density 
field and the matter power spectrum is very small. The small corrections we have found are 
well fit away from A; = by a polynomial. We list best fit polynomial parameters in Tabled 
in Appendix B. The amplitude of this correction in our fiducial cosmology should provide a 
conservative estimate of the uncertainty in the c orrection as the c osmo logy is varied within 
the space allowed by the latest WMAP analysis iKomatsu et al.l (120081 ): ~ 1% for k < 0.1 
and ~ 4% between A; = 0.1 and k = 0.2. Finally, to recover the linear power spectrum, 
one must know the amplitude of the BAO suppression and the degree of nonlinearity in the 
matter power spectrum as a function of the cosrnological parameter s. Several groups are 



addressing this issue (e.g., lEisenstein et al.ll2007bl : iHabib et al.ll2007l ). 



6.2.2. Covariance Matrices 

We use the scatter in the power spectrum measurements in the simulations to estimate 
the covariance matrix for our mock catalogs. We fit the measurements of the four-point 
function to a physically motivated model: 

B,, = AR^P{h)P{k,)6lc{z) (31) 
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Fig. 7. — Upper left. The ratio of the monopole (i.e., angle-averaged) LRG P{k) to the real 
space nonhnear dark matter spectrum for several MID subsamp les: central LRGs only (solid 
curves); all LRGs (dotted curves); all LRGs after undergoing iTegmark et al.l (l2006l ) FOG 
compression (dash-dot curve); and CiC groups (dashed curve). For the first two subsamples 
we also show the real space P{k), which have lower amplitude as — > 0. Upper right: 
The real space P{k) for the central and satellite LRGs minus b^^i times Pcenik), the power 
spectrum of the central LRGs in real space, for hrei = 1-037, 1.042, 1.047. Lower left: The 
solid curve shows the ratio of the central LRGs P{k) in redshift space to the central LRGs 
P{k) in real space, while the dashed curve shows the ratio for Pcic{k)- Lower right: The 
redshift space monopole P{k) for the central and satellite LRGs minus 6^^; times Pcen{k), the 
power spectrum of the central LRGs in redshift space, for brei = 1-037, 1-042, 1.047. Error 
bars are similar for all curves in a figure, and so we show only one set for clarity. 
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Fig. 8. — Same as Fig. [7] upper left panel, but for the NEAR sample. 
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Fig. 9. — Same as Fig. [7] upper left panel, but for the FAR sample. 
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Fig. 10. — We separately fit o ur NEAR, MID , and FAR FOG-compressed samples to the 
nonlinear correction model of iTegmark et al.l (l2006l ) in Eqn. [6l Best fit parameters are 
^NEAR = 5.61, Qnl,near = 17.4 (dotted curve); = 6.15 , Qnl,mid = 18.4 (dashed 

curve); and bp^j^ = 6.67, Qnl,far = 23.3 (sohd curve). Each subsample is scaled by the 
input power spectrum and the best fit 6^ value to accentuate the variation in the power 
spectrum shapes for these three samples. The NEAR and FAR data points and model 
cu rves have been s hifted by -/+ 0.0013 h Mpc~^ for clarity. The best fit Qnl vanilla model 



m 



Tegmark et al.l (120061 ). A = 1.4 and Qnl = 31 in Eqn. HI is shown by the long dashed 



curve. We have adjusted the normalization to agree at k ^ 0.05. The do t-dashed curve 
shows the ratio of the best fit linear power spectrum of iTegmark et al.l (120061 ) and the linear 
power spectrum adopted in this work, demonstrating that Qnl is degenerate with changes 
in the spectral shape, so that cosmological parameter information is lost to Qnl- 
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C,,,tot = ^ {p{h)P{kj) + '^y 5^^+(3totB,, (32) 

= -j- (P\h) + "'•^'^^^^^^ + — ^ + (33) 

A^i V n^NmocksJ 

- iv^l n + n^(iV_,.-i)J^-+^-^-^- (2^^ 

Note we use the measured P{ki) rather than the hnear or model P{ki) in Eqns. [MllMl In the 
reduced covariance we expect the Poisson noise to be suppressed by l/Nmocks, the number of 
mock catalogs we produce for each TPM simulation. Eqn. [T9] implies that anoD+o^red = '^cttot 
for atot ~ 1 and 13 hod + [3red = Ptot- Best fit values for each mock catalog type and each 
redshift subsample are listed in Table HI where the fit included all modes with k < 0.4, and 
in Table [5], where the fit included all modes between k = 0.05 and k = 0.2. We present both 
in order to assess whether the non-Gaussian terms are growing with k. In Fig. [11] we show 
the measured and model Cij values along the diagonal as well as three rows in the matrix for 
the MID reconstructed halo density field covariance matrix; the agreement with the model 
for other rows and redshift subsamples are similar. 

We find atot ~ 1 for all mocks and redshift subsamples, indicating that the standard shot 
noise contribution to the covariance is approximately correct for our mock catalogs. However, 
the relative distribution of this cross term between Cij^noD and Cij^red varies; a hod is large 
for the FAR sample, where the number density is ~ half that in the MID and NEAR sample, 
and is also larger in the samples including satellites relative to the samples with only central 
objects. All three matrices have significant off-diagonal terms. f3tot is larger than unity for 
all samples, which may be expected since the LRGs are biased tracers, so the DC mode 
variance will increase by the factor b^i^c- However, the amplitude does not scale with 6^, but 
is much larger for the FAR sample. The best fit parameters do not change between real and 
redshift space for the same sample, but Ptot is larger for the samples including satellites. The 
reconstructed halo density field best fit parameters are consistent with the parameters of the 
sample of central galaxies, and the FOG-compressed sample is consistent with the sample 
of central and satellite galaxies. Ptot decreases slightly when fitting only to the k = 0.05 to 
k = 0.2 results, but the effect is small and so we adopt the parameters reported in Table [5] 
to carry out our covariance matrix calculations for the error bars on PLRcik)- 

There is a modest decrease in the normalized diagonal covariance Cii/P{kiY (~ 5%) 
when satellite galaxies are included. However, this decrease is offset by the fact that the 
P^(fcj) has a larger nonlinear component, so not all of the available information will be about 
the linear component; furthermore, the large nonlinear correction introduces uncertainty in 
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extracting the linear component. Moreover, l3tot is larger for samples which include satellite 
galaxies, so the bandpowers will be more highly correlated for those mocks. We therefore 
conclude that the reduction in large scale bias by eliminating the satellite galaxies and recon- 
structing the halo density field is offset in the error budget by smaller off-diagonal covariance 
and smaller uncertainties in the nonlinear correction. 



6.2.3. The Redshift Space Quadrupole 



In linear theory the redshift space power spectrum can be decomposed into a monopole 
quadrupole, and hexadecapole: 

1 dlnD 



Pjk) 



bgai d\na 
P{k) 



(35) 
(36) 



where Lj is the Legendre polynomial of order i , and fi^ = k ■ s, with s the direction along 
the LOS. The method of iTegmark et al.l (120061 ) assumes Eqn. [36] with /3 independent of k 
to obtain the best estimate of the real space power spectrum from Pgg, Pg^, and P^. In 
Fig. [12] we examine the structure of the redshift space distortions in our mock catalogs. 
The first three panels show the quadrupole to monopole ratio for central LRGs (solid), the 
reconstructed halo density field (dashed), and the FOG-compressed density field (dash-dot). 
These are all similar and show a modest increase in the qua drupole to monopole r atio with 
k. This demonstrates that the FOG compression technique of iTegmark et al.l (l2006l ) success- 
fully removes the effects of the FOGs induced by satellite galaxies. The oscillations probabl y 
result from the extra suppression of BAO features in redshift space (lEisenstein et al.ll2007bl ). 
The central and satellite LRG sample (dotted) quadrupole to monopole ratio lies below the 
expected value at all k values, indicating that the FOG suppression of power is evident even 
in the linear regime. Therefore a sample where the effects of satellite FOGs have not been 
removed by either halo density field reconstruction or FOG compression will provide a biased 
estimate of j3. 



The bottom right panel shows the ratio of the redshift space monopole to the real space 
power spectrum for the sample of central galaxies only (solid) and central and satellite galax- 
ies (dotted). The suppression of power by the satellite FOGs is also evident here. For central 
LRGs only, the ratio is consistent with the expected value from Eqn. [36], and between k = 
and k = 0.2 this ratio varies by only ~ 3%. While there is a clear scale dependence in the 
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kj h/Mpc 



Fig. 11. — We plot the normalized covariance matrix elements Cij/ P{ki)P{kj) scaled by 
a/ NiNj, the number modes in bins i and j. We plot our measurements of Cij^tot (dotted 
curve), Cij^red (dashed curve), and Cij^HOD (dot-dashed curve). The solid lines show the 
model fits from Table HI The top left panel shows the diagonal elements, and the other 
panels show cross sections of the covariance matrix with fcj = 0.09, 0.15, and 0.3 h Mpc~^. 
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quadrupole to monopole ratio for the central LRG sample, the large bias of LRGs means that 
the monopole redshift space power spectrum is nearly insensitive to the redshift space non- 
linearities in the halo density field. The upper dotted curve shows Ps,sat/ Preai,cen x ^cen/^'saf 
The increase in power in real space when the satellites are included is larger than the sup- 
pression of power in redshift space by their FOGs, so that the satellite monopole spectrum 
has more power at high k than the real space central LRGs. 



7. Conclusions 



This paper introduced and tested our algorithm for using the halo density field to es- 
timate the underlying matter power spectrum. We found the nonlinear correction between 
this field and the underlying matter density field to be both smaller and more robust to vari- 
ations in the e ffects of satellite galax ies than both the FOG-compressed density field used in 
the analysis of iTegmark et al. and the redshift space monopole power spectrum used 

in the analysis of iPercival et al.l (120071 ). 



The parameters of our simulation set were selected to provide accurate two point halo 
density and velocity statistics to assure accurate representations of FOG features in our 
mock catalogs. To our knowledge, this is the first attempt to study the detailed effects of 
FOG treatment on the LRG power spectrum. We find that the FOG treatment can affect 
bandpowers even in the linear regime, k < 0.1. We first examined the nonlinear matter 
power spectrum of our 42 simulations. The ratio of the nonlinea r to input matter power 
spectra is well described by a smearing of the BAOs as modeled by lEisenstein et al.l (l2007bl ) 
and a smooth increase in power with k that can be fit by a third order polynomial out to 
k = 0.4 h Mpc^^ or second order polynomial out to k = 0.2 h M pc~^. We detect a substan- 
tial deviation from the predictions of halofit (jSmith et al.ll2003l ). 



Fig. [7] demonstrates the main point of this work: satellite galaxies systematically alter 
the shape of the power spectrum even at k < 0.1 h Mpc^^. Extraction of cosmological infor- 
mation from the broad band shape of the power spectrum is already limited by systematics 
( ISanchez fc Colell2007l ) which we (and others) suggest can be attributed primarily to differ- 
ences in the satellite contribution to the power spectrum. In t his paper we demonstrate that 
while the FOG compression scheme in iTegmark et al.l (120061 ) exacerbates these issues and 
requires a large nonlinear correction, the FOG features in the density field can be used to 
reconstruct the halo density field with high fidelity. The power spectrum of this field deviates 
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k h/Mpc 



Fig. 12. — The bottom right panel shows the ratio of the redshift space monopole to the 
real space power spectrum for the MID sample. The straight lines show the prediction of 
Eqn. [36] for both samples. We fit the bias in Eqn. [35] to the large scale ratio PlrgI Pdm- 
The solid curve is for the sample of central galaxies only {been = 2.11) and the dotted curve 
with error bars includes both central and satellite galaxies {bgat = 2.23). The upper dotted 
curve shows the ratio of the redshift space monopole of the central and satellite galaxies to 
the real space power spectrum of the central objects (scaled by fccen/^Lt)- The other three 
panels show the quadrupole to monopole ratio in redshift space. Error bars are only shown 
for the central galaxy sample (solid curve) for clarity. The central and satellite sample is 
shown by the dotted curve, the reconstructed halo density field as the dashed curve, and 
the FOG-compressed sample as the dash-dot curve. The expected values from Eqn. [36] are 
shown as the straight lines for the central (solid) and central and satellite (dotted) samples. 
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from the dark matter power spectrum at the < 4% level for k < 0.2 h Mpc~^ and < 1% level 
for k < 0.1 h Mpc~^, where cosmological analyses usually restrict themselves. Moreover, we 
have shown that this correction changes only slightly between the NEAR, MID, and FAR 
LRG reconstructed halo density fields, while the FOG compressed mocks have a much larger 
variation between samples. Therefore, we can hope to push cosmological analyses to larger 
k using the reconstructed halo density field as a tracer of the underlying matter density field 
fluctuations — particularly for galaxy samples like the LRGs which are spread over a large 
redshift range and are not volume-limited, and thus have substantial variation in the satellite 
contribution to the power spectrum with redshift. 



While we have not addressed the variation of the nonlinear correction to the recon- 
structed halo density field as a function of cosmological parameters, we have designed the 
form of our corre ction to minimize the variation with cosmology. Other researchers (e.g.. 



Habib et al.ll2007l ) are studying the dark matter power spectrum as a function of cosmology. 



We expect that our nonlinear correction Plrg/ Pdm will remain small (of order 4% below 
k = 0.2 h Mpc~^) as the cosmology is varied, so the variation of this small correction should 
be even smaller. Therefore, instead of introducing a nuisance parameter for the nonlinear 
correction, we propose that the amplitude of the correction should be taken as the error on 
its value in cosmological parameter analyses, or relatively strong priors on the amplitude of 
the correction to the reconstructed halo density field be introduced. 



In this work we have also investigated the properties of the power spectrum covariance 
matrix for the dark matter, as well as the mock galaxy catalogs divided into 6 different 
samples: central galaxies in real and redshift space, central and satellite gala xies in real and 



redshi ft space, our reconstructed halo density field in redshift space, and the iTegmark et al. 



(120061 ) FOG compressed galaxy density field in redshift space. All of these samples were 
well-modeled by a diagonal matrix with the usual Gaussian a nd Poisson shot noise terms 
plus the beat-coupling term presented in [Hamilton et al.l (120061 ). We expect that this will be 
a useful model for fitting the survey covariance matrix, where P{kh) is replaced by Piksurvey) 



with k 



survey 



determined by some effective survey size. 



Finally we examined the structure of the redshift space distortions as a function of k 
using the quadrupole. Both the reconstructed halo density field and the FOG-compressed 
mock LRG samples reproduce the modest k dependence of the halo density field quadrupole 
to monopole ratio. Since the LRGs are so highly biased, this scale dependence causes a < 3% 
deviation in the redshift space monopole to real space power spectrum ratio out to A; = 0.2. 
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When satellites are included in the sample without FOG compression, the quadrupole to 
monopole ratio is lower than the linear value for the entire k range accessible in our simula- 
tions and will significantly bias the estimate of f3. 
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9. Appendix A: Resolution Study Results 



The goal of the resolution study detailed in Chapter 3 of iReidl ( 120081 ) was to determine 
the minimum iV-body simulation mass resolution satisfying two criteria. 



The minimum halo mass is small enough to accommodate broad Ncen{M) functions 
(i.e., large values of aiogu)- The width of this function is dictated by the observed large- 
scale clustering amplitude, either in PiRGik) or Wp{rp), which determines asbiRG- 

The finite mass resolution does not systematically bias the small scale n-point halo 
statistics that must be acc urate in or der to test our method of halo density field 
reconstruction. We argue in lReidI (120081 ) that agreement of several two-point clustering 
and velocity statistics and the CiC group multiplicity functions is sufficient for accurate 
mock FOG structures. Furthermore, the small observed deviations in the halo mass 
function at low M can be absorbed by the parameters of Ncen{M). 



For the cosmological parameters of our mock catalogs and the observed LRG clustering 
strength, the first condition places the more stringent requirement on the simulation mass 
resolution. To establish the level of systematics relevant to the second condition, we com- 
pare the results of one iV-body simulation with our adopted parameters Lhox = 558 h^^ 
Mpc and Np^med = 512'^ {Mpmed = 1-43 x IO^^Mq) with a higher resolution simulation 



p,high 



640^ M 



p,high 



7.33 X 10^°Mq) and a lower resolution simulation {Nj 



p,low 



3843, 



Mp^iow = 3.39 X IO^^Mq) of the same initial conditions. We highlight the findings of this 
comparis on below for the simulation outputs near our zmid = 0.342, and refer the interested 
reader to iReidI (120081 ) for a more thorough discussion. 



This preprint was prepared with the AAS macros v5.2. 
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9.1. Resolution Effects in the Halo Mass Function and the Halo Occupation 

Distribution 

The mass functions from the low, medium, and high resolutio n simulations are in excel- 



lent agreement with one another and within a few percent of the [Tinker et al.l (120081 ) mass 
function. There is a slight increase (~ 5%) in the halo mass function as the aggressive mass 
limit 50Mp of our SO halo catalogs is reached; this should be kept in mind in future analyses, 
but is not inherently problematic since slight changes in the mass function in this regime are 
degenerate with small variations in the form of Ncen{M). 



We produce LRG mock catalogs at aiogM = 0.6, the value that matches the observed 
LRG clustering amplitude, in the low, medium, and high resolution simulations. For the high 
resolution simulation, only 2.6% of the mock LRGs occupy halos below the halo catalog mass 
limit of the medium resolution simulation, and only 11% are below lOOMp^med- Therefore, 
the abrupt cutoff in Ncen{M) induced by our halo catalog mass limit will show minimal dif- 
ferences with a mock catalog derived from a halo catalog extending to a substantially lower 
mass limit. In addition, even if there are slight systematics in the pair statistics of a halo cat- 
alog in the ~ 50 — lOOMp^med mass range, they will be suppressed by at least a factor of ~ 10. 



9.2. Two-point clustering and velocity statistics 

There is a slight but systematic decrease in power with k as the simulation resolution 
decreases. This is a 0.2% effect at A; = 0.2 h Mpc~^ between the medium and high resolution 
simulations, well within the statistical errors and modeling uncertainties between the LRGs 
and dark matter density fields. 

We split the halos into two mass bins: 50 — 99Mp^rned and M > lOOMp^med- This corre- 
sponds to 98 — 195 particles per halo for the low mass bin of the high resolution simulation. 
We did not detect any significant resolution-dependent differences in P{k) and C,{r), other 
than a < 1% change in the large scale bias of the halo samples with respect to the mat- 
ter. Fig. [13] shows that the difference between the medium and high resolution simulation 
^jj(r), for both auto- and cross-correlations of the low and high mass halo bins, is consis- 
tent with at small scales {r ^ 15 Mpc); this is true out to 150 Mpc. Fig. 

shows the velocity correlation statistic (vj(r) ■ Vj(r')) / ^ (vf ) ^v|) for halo mass bins i and 
j. We find < 1% level agreement between the medium and high resolution simulations. 
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In Fig. [H] we also examine the relative halo velocities along the vector separating the ha- 
los: ((vj(r) — Vj(r')) • (r — r')/|r — r'|) and find no evidence for systematic bias between the 
medium and high resolution simulations. 



9.3. Mock catalog statistics: two point clustering statistics and multiplicity 

bias b{nsat) 



In iReidI (l2008l ) we compare both the power spectrum P{k) and the projected correlation 
function Wp{rp) of the medium and high resolution simulations, for which we find systematic 
deviations smaller than ~ 1/6 the statistical error induced by the stochasticity of the halo 
occupation. We quantify the difference between the true one-halo group density field and the 
CiC group density field by the multiphcity bias bijisat) = Ncicinsat) I Ntmeinsat) ■ We show 
that b{nsat) systematically varies by less than the amount of variation induced by the halo 
occupation stochasticity for mock catalogs with aiogM = 0.8, which pushes to lower mass 
limits than our better matched aiogM = 0.6 catalogs. 



We conclude that the selected mass resolution of our simulations is well-optimized to 
produce a large simulation volume over several realizations and with well-controlled system- 
atic biases induced by the finite mass resolution. 



10. Appendix B: Power Spectra and Covariance Matrix Best Fit Parameters 
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Fig. 13. — Upper left: r^ij{r) on small scales for the low mass halos auto-correlation (solid), 
high mass halos auto-correlation (dashed), and their cross correlation (dotted). Also shown 
are {ihigh{i^) — imed{r)) / C.highi'f^) Computed in bins of size Ar = 0.5 Mpc for the low mass 
halos auto-correlation (upper right) , high mass halos auto-correlation (lower right) , and their 
cross correlation (lower left). The agreement for each correlation extends to r = 150 
Mpc (not shown). 
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Fig. 14. — Upper left: The normalized halo-halo velocity correlation function, 
(vj(r) ■ Vj(r')) / -y/(v^)~(v|y, vs |r — r'| for the low mass halos auto-correlation (solid), high 
mass halos auto-correlation (dashed), and their cross correlation (dotted) binned in in- 
tervals of 0.5 Mpc in |r — r'|. Upper right: The fractional difference in the normal- 
ized halo-halo velocity correlation between the high and medium resolution simulations 
binned in intervals of 5 Mpc in |r — r'|. Lower panels: Same as the upper panels, 
but for the halo-halo relative velocity correlation function along the halo separation vector, 
((v,(r)-v,(r'))-(r-r')/|r-r'|). 
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Table 3. Fits to nonlinear power spectra, k^axjit and ksAo are in units of h Mpc~^. For 

PDhiik) / Pic{k) we fit to Eqn. [301 When kmaxju = 0.2 we hold = 0. Some of the 
redshift space mock catalog fits are not well-behaved at very small k's, and the corrections 
to a constant should be suppressed at < 0.05. This also makes deviate from the large 
scale bias, given in the last column and fit using only k < 0.06. Psat,z{k) denotes the power 
spectrum of a catalog containing central and satellite galaxies in redshift space, Pcic,z{k) 
denotes the power spectrum of the reconstructed halo de nsity field in redshift space, and 
Pteg,zik) denotes an FOG compressed following iTegmark et al.l (120061 ). 



sample 


function 


^inax ,fit 


ksAO 


10 


ai 


(12 


0.3 


"lrg 


NEAR 


PDM{k)/Pic(k) 


0.2 


0.12 


1.025 


-1.109 


9.296 






NEAR 


PDM{k)IPic(k) 


0.4 


0.13 


1.027 


-1.305 


12.05 


-9.136 




NEAR 


PsatAk)/PDM(k) 


0.2 




5.14 


5.38 


-15.96 




5.38 


NEAR 


PsatAk)/PDM(k) 


0.4 




5.01 


7.91 


-28.72 


15.08 


5.38 


NEAR 


PciC,Ak)/PDM{k) 


0.2 




4.74 


2.11 


-12.80 




4.83 


NEAR 


PciC,z{k)/PDM{k) 


0.4 




4.66 


3.68 


-20.45 


8.40 


4.83 


NEAR 


PtegAk)/PDM{k) 


0.2 




5.15 


4.85 


7.78 




5.42 


NEAR 


PtegAk)/PDM{k) 


0.4 




5.03 


6.44 


6.41 


-20.33 


5.42 


MID 


PDM{k)IPic(k) 


0.2 


0.13 


1.025 


-1.042 


8.614 






MID 


PDM{k)/Pic{k) 


0.4 


0.14 


1.026 


-1.199 


11.06 


-8.426 




MID 


PsatAk)/PDM(k) 


0.2 




5.67 


6.89 


-17.9 




5.99 


MID 


PsatAk)/PDM(k) 


0.4 




5.54 


9.21 


-28.3 


8.54 


5.99 


MID 


Pcic(k)/PDM(k) 


0.2 




5.30 


3.43 


-14.7 




5.46 


MID 


Pc,c{k)/PDM{k) 


0.4 




5.21 


5.13 


-21.8 


4.35 


5.46 


MID 


Pteg{k)/PDM{k) 


0.2 




5.68 


6.11 


8.23 




6.03 


MID 


Pteg{k)/PDM(.k) 


0.4 




5.58 


7.27 


10.4 


-28.4 


6.03 


FAR 


PDM{k)/Pic{k) 


0.2 


0.13 


1.025 


-0.994 


8.144 






FAR 


PDM{k)/Pic{k) 


0.4 


0.15 


1.025 


-1.126 


10.39 


-7.945 




FAR 


PsatAk)/PDM(k) 


0.2 




6.12 


9.72 


-19.26 




6.58 


FAR 


PsatAk)/PDM(k) 


0.4 




5.89 


13.96 


-37.46 


11.39 


6.58 


FAR 


Pcic(k)/PDM(k) 


0.2 




5.79 


4.36 


-13.97 




6.00 


FAR 


Pcic(k)/PDM{k) 


0.4 




5.64 


6.94 


-23.48 




6.00 


FAR 


Pteg(k)/PDM(k) 


0.2 




6.16 


7.90 


26.43 




6.63 


FAR 


Pteg(k)/PDM(k) 


0.4 




5.96 


10.49 


26.82 


-45.03 


6.63 
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Table 4. Fits for the mock catalog covariance matrices using Eqns. |3T]- [311 ^ill 
bandpowers between k = 0.01 and k = 0.4 are included in the fits, 'cen, real/redshift' 
denotes a sample of only central LRGs in real/redshift space; 'sat, real/redshift' denotes a 
sample with both central and satellite galaxies in real/redshift space; 'CiC, redshift' 
denotes the reconstructed halo density field using the Ci C method, and ' Teg, r edshift' 
denotes the FOG-compressed sample according to iTegmark et al.l (120061 ). 



Catalog 


ohod 


^red 


2atot 


Phod 


0red 


Ptot 


NEAR cen, real 


1 


07 


1.10 


2 


00 


0.26 


0.94 


1.26 


NEAR sat, real 


1 


30 


0.89 


1 


98 


0.36 


1.58 


1.97 


NEAR cen, redshift 


1 


11 


0.99 


2 


02 


0.26 


0.95 


1.29 


NEAR sat, redshift 


1 


27 


0.63 


2 


01 


0.29 


1.28 


1.64 


NEAR CiC, redshift 


1 


11 


0.92 


2 


02 


0.26 


0.94 


1.30 


NEAR Teg, redshift 


1 


39 


0.74 


1 


98 


0.35 


1.60 


1.99 


MID cen, real 


1 


14 


1.0 


2 


00 


0.34 


1.21 


1.62 


MID sat, real 


1 


32 


0.77 


1 


99 


0.42 


1.71 


2.20 


MID cen, redshift 


1 


17 


0.93 


2 


01 


0.34 


1.19 


1.63 


MID sat, redshift 


1 


30 


0.61 


2 


00 


0.36 


1.47 


1.92 


MID CiC, redshift 


1 


16 


0.86 


2 


01 


0.34 


1.18 


1.63 


MID Teg, redshift 


1 


41 


0.58 


1 


97 


0.41 


1.73 


2.20 


EAR cen, real 


1 


65 


0.50 


2 


05 


1.32 


1.62 


3.13 


FAR sat, real 


1 


85 


0.38 


2 


07 


1.71 


2.67 


4.50 


FAR cen, redshift 


1 


66 


0.42 


2 


05 


1.25 


1.64 


3.13 


FAR sat, redshift 


1 


79 


0.16 


2 


07 


1.37 


2.08 


3.70 


FAR CiC, redshift 


1 


66 


0.37 


2 


06 


1.25 


1.69 


3.21 


FAR Teg, redshift 


1 


90 


0.13 


2 


06 


1.54 


2.53 


4.19 
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Table 5. Fits for the mock catalog covariance matrices using Eqns. |3T]- [341 all 
bandpowers between k = 0.05 and k = 0.2 are included in the fits. Samples are the same as 

in Table H 



Catalog CtHOD Qred 'iatot PhOD Pred Plot 



NEAR cen, real 


1, 


,05 


0, 


,85 


1, 


,96 


0, 


,23 


0, 


,86 


1, 


,14 


NEAR sat, real 


1, 


.22 


0, 


,62 


1, 


,94 


0, 


,29 


1, 


,43 


1, 


.76 


NEAR cen, redshift 


1, 


,07 


0, 


,71 


1, 


,90 


0, 


,23 


0, 


,94 


1, 


.23 


NEAR sat, redshift 


1, 


,22 


0, 


,51 


1, 


,86 


0, 


,26 


1, 


,33 


1, 


.66 


NEAR CiC, redshift 


1, 


,07 


0, 


,68 


1, 


,90 


0, 


,24 


0, 


,97 


1, 


.27 


NEAR Teg, redshift 


1, 


,28 


0, 


,45 


1, 


,86 


0, 


,29 


1, 


,47 


1, 


.82 


MID cen, real 


1, 


,08 


0, 


,84 


1, 


,96 


0, 


,28 


1, 


,14 


1, 


.49 


MID sat, real 


1, 


,21 


0, 


,62 


1, 


,93 


0, 


,35 


1, 


,60 


2, 


.02 


MID cen, redshift 


1, 


,11 


0, 


,74 


1, 


,93 


0, 


,30 


1, 


,20 


1, 


.58 


MID sat, redshift 


1, 


,23 


0, 


,54 


1, 


,88 


0, 


,33 


1, 


,55 


1, 


.96 


MID CiC, redshift 


1, 


,10 


0, 


,72 


1, 


,91 


0, 


,29 


1, 


,24 


1, 


.61 


MID Teg, redshift 


1, 


,29 


0, 


,46 


1, 


,88 


0, 


,36 


1, 


,66 


2, 


.10 


EAR cen, real 


1, 


,57 


0, 


,36 


2, 


,00 


1, 


,08 


1, 


,66 


2, 


.87 


FAR sat, real 


1, 


,68 


0, 


,21 


2, 


,00 


1, 


,35 


2, 


,51 


3, 


.98 


FAR cen, redshift 


1, 


,57 


0, 


,26 


1, 


,99 


1, 


,06 


1, 


,70 


2, 


.94 


FAR sat, redshift 


1, 


,67 


0, 


,13 


1, 


,98 


1, 


,20 


2, 


,28 


3, 


.65 


FAR CiC, redshift 


1, 


,58 


0, 


,26 


2, 


,00 


1, 


,06 


1, 


,84 


3, 


.09 


FAR Teg, redshift 


1, 


,73 


0, 


,06 


1, 


,99 


1, 


,31 


2, 


,47 


3, 


.93 



